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Abstract 

This paper deals with fractional-order controllers. 
We outline mathematical description of fractional 
controllers and methods of their synthesis and ap- 
plication. Synthesis method is a modified root locus 
method for fractional-order systems and fractional- 
order controllers. In the next section we describe how 
to apply the fractional controller on control systems. 
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1 Introduction 

The real objects are generally fractional how- 
ever, for many of them the fractionality is very low. 
A typical example of a non-integer (fractional) or- 
der system is the voltage-current relation of a semi- 
infinite lossy RC line or diffusion of the heat into a 
semi-infinite solid, where heat flow is equal to the 
half-derivative of temperature p6[ : 



dP''T{t) 



lit). 



Their integer-order description can cause significant 
differences between mathematical model and the real 
system. The main reason for using integer-order 
models was the absence of solution methods for 
fractional-order differential equations. Recently, im- 
portant achievements [0, |, |, |, ^ were 
obtained, which enable to take into account real or- 
der of dynamic systems. 

PID controllers belong to dominating industrial 
controllers and therefore are objects of steady ef- 
fort for improvements of their quality and robustness. 
One of the possibilities to improve PID controllers 
is to use fractional-order controllers with non-integer 
derivation and integration parts. For fractional-order 
systems the fractional controller CRONE ||] has 
been developed, PD^ controller ||l and the PI^D^ 
controller ||l2|, |l^ has been suggested. 
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Before the actual design of controllers for dynam- 
ical systems it is necessary to identify these sys- 
and then determine their dynamical 
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tems 

properties. The dynamical properties of the system 
under observation can be expressed mathematically 
and graphically in the form of various characteristics. 
These characteristics can be determined from the dif- 
ferential equation or transfer function via computa- 
tion or experimentally by exciting the system from 
the equilibrium. Computation of the transfer char- 
acteristics of the fractional-order dynamical systems 
has been the subject of several publications, e.g. by 
numerical methods Q , as well as analytical methods 
[ p^ . In the synthesis of a controller its parameters 
are determined according to the given requirements. 
These requirements are, for example, the stability 
measure, the accuracy of the regulation process, dy- 
namical properties etc. The check whether the re- 
quirements have been met can be done with a sim- 
ulation on the control circuit model. There are a 
large number of methods for the design of integer- 
order controllers, but the situation is worse in the 
case of fractional-order controllers where the meth- 
ods are only being worked out. One of the methods 
being developed is the method (modification of roots 
locus method) of dominant roots jlO|, based on the 
given stability measure and the damping measure of 
the control circuit. 



2 Basic mathematical tools for 
fractional calculus 

The fractional calculus is a generalization of integra- 
tion and derivation to non-integer order operators. 
The idea of fractional calculus has been known since 
the development of the regular calculus, with the first 
reference probably being associated with Leibniz and 
L'Hospital in 1695. At first, we generalize the differ- 
ential and integral operators into one fundamental 
operator Z?" which is known as fractional calculus: 
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^{a) > 0, 
^(a) = 0, 
m^a) < 0. 



The two definitions used for the general frac- 
tional differintegral are the Griinwald definition and 



the Riemann-Liouville (RL) definition 
Griinwald definition is given liere 
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where [x] means the integer part of x. The RL defi- 
nition is given as 



1 



-dr, (2) 



for (n—l<a<n) and where T{x) is the well known 
Euler's gamma function. 

The Laplace transform method is used for solving 
engineering problems. The formula for the Laplace 
transform of the RL fractional derivative has the 
form [Pi: 



(3) 
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for (n — 1 < a < n). 

For numerical calculation of fractional-order 
derivation we can use the relation (^) derived from 
the Griinwald definition (|l|). This relation has the 
following form: 

Nit) 

where L is the "memory length", h is the step size of 
the calculation, 



Nit) 



[x] is the integer part of x and bj is the binomial 
coefficient: 



bo = 1, 
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For the solution of the fractional-order differential 
equations (FODE) most effective and easy analytic 
methods were developed based on the formula of the 
Laplace transform method of the Mittag-Leffier func- 
tion in two parameters A two-parameter func- 
tion of the Mittag-Leffier type is defined by the series 
expansion: 



k=0 
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(a,/?>0). (6) 



The Mittag-Leffier function is a generalization of ex- 
ponential function and the exponential function is 
a particular case of the Mittag-Leffier function. Here 
is the relationship given in |13||: 
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3 Fractional - order closed 
control loop 

We will be studying feed-back control system with 
unit gain in the feed-back loop (Fig.l), where Gr{p) 
is the controller transfer function, Gg (p) is controlled 
system transfer function, W{p) is an input, E{p) is 
an error, U{p) is output from controller and Y{p) is 
output from system. The transfer function of close 
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Fig.l: Feed - back control circuit 
feed-back control circuit (Fig.l) has the form: 

_ Gr{p)Gs{p) 

"^'^P'-Wip) l + Grip)Gsip)- 

3.1 Fractional - order controlled 
systems 

The fractional-order controlled system will be repre- 
sented with fractional model with transfer function 
given by the following expression [Q: 



Gsip) 
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where f3k, (fc = 0, 1, 2, . . . , n) are generally real num- 
bers, (3n > Pn-1 > • ■ • > /9i > /So > and Ofc 
(k — 0,1, . . . ,n) are arbitrary constants. In the time 
domain, the transfer function (||) corresponds to the 
n-term FODE with constant coefficients 

a„ D^"yit) + . . . + ai D^'y{t) + ao D^"y{t) = u{t). 

(9) 

Identification methods |l^ for determination of 
the coefficients ak and (3k, (k = 0,1, . . . ,n) were de- 
veloped, based on minimization of the difference be- 
tween the calculated [y'^) and experimentally mea- 
sured {y'^) values 



Q = 
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where M is the number of measured values. 

For the analytical solution of the n-term FODE 
we can write formula in general form p!3]: 
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where E\,^{z) is the Mittag-Leffler function in two 
parameters (o), 



{j + n)\ 



for (n = 0, l,2,...)• 
For the numerical solution of the n-term FODE 
we can write formula in general form: 



for A: = 1, 2, 3, . . ., y(0) = 0, = 0, where u{k) is a 
function on the right side of the differential equation. 

3.2 Fractional - order controllers 

The fractional-order controller will be represented by 
fractional-order PI^D^ controller with transfer func- 
tion given by the following expression p4|: 



(10) 



where A and 5 are an arbitrary real numbers (A, 5 > 
0), if is amplification (gain), Ti is integration con- 
stant and Td is differentiation constant. In the time 
domain equation (^) has the form: 

u{t) = Ke{t) + TiD-^e{t) + TdDfe{t). (11) 

Taking A = 1 and S = 1, we obtain a classical 
PID controller. If A = (T^ = 0) we obtain a PD^ 
controller, etc. All these types of controllers are par- 
ticular cases of the PI^D^ controller. The PI^D^ 
controller ( pj] ) is more flexible and gives an opportu- 
nity to better adjust the dynamical properties of the 
fractional-order control system. 

4 Synthesis of fractional - 
order controllers 

For the design of fractional-order controller a new 
method was suggested based on the dominant roots 
principle. This method is based on from poles distri- 
bution of the characteristic equation in the complex 
plane (Fig.2). 

Values of dominant roots are designed for the qual- 
ity requirement of the control circuit. Their signifi- 
cance is, that dominant roots are defined for stability 
measure St and damping measure Ti. Under there 
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Fig.2: Roots in the complex plane 



conditions the complex conjugate roots satisfy the 
equation: 

Pi,2 — -r ± iuj. (12) 

The parameters of controller were set up so that 
other poles were from dominant remote to the left 
side. The parameters design of the fractional-order 
controller can be divided into two stages: 

1. Design of parameter K 

Proportional parameter K influence the value 
of static deviation Et [%] , control time T^. [s] , 
and overshoot Pr [%]■ Generally, with increased 
value of K, control time Tr [s] is decreasing and 
static deviation Et [%] is lower: 

K > [lOO/Et] - ao. 

2. Design of parameters Td, S, Ti, A 

We define the required stability measure St = r 
and damping measure Ti = ^. This requirement 
is satisfied by the complex conjugate roots (^2|). 
We will similarly use characteristic equation as 
the classical root locus method [|| . The charac- 
teristic equation of fractional-order control loop 
(0) has the form: 



Gr{p)Gs{p) + 1 = 0. 



(13) 



After substitution of the fractional-order con- 
troller transfer function (|l^) and the fractional- 
order controlled system transfer function (^) and 
after some manipulation wc obtain characteris- 
tic equation (O) in the following form: 



Y.akP^" + iK + Tp-^ + TdP^) = 0. (14) 

k=0 



This algebraic equation is in the general form 
for the n-term fractional-order controlled system 
and fractional-order PI^D^ controller. Solution 
of this equation for known poles and parameter 
K gives unknown parameters Tcj, (5, T^, A. The 
algebraic equation is solved in the space 
of complex variable. Result are parameters of 
fractional-order controller for required stability 
measure and damping measure. 

5 Control algorithm for 

fractional - order controllers 

The control algorithm was designed according to the 
control scheme in Fig. 1 . The position algorithm uses 
discrete time steps k and consists of the following 
steps 0: 

1. Filtering of required value: 

w*{k) =w*{k-l) + 0.5{w{k) - w*{k - 1)), 
where w{k) is required value. 

2. Calculating the control error: 

e(fc) = w*{k) - y{k), 

for discrete time step (fc = 1, 2, . . .), where e(fc) 
is regulation error and y{k) is measured value. 

3. Determination of control value: 



uik) = K e{k) + ^Y1 - J') + 

j=v 

for discrete time step (fc = l,2,...), where T is 
the length of time step (sample period). The 
binomial coefhcients dj and q,- were calculated 
from the recurrent equation (j^. The numerical 
algorithm (^) requires to store the whole history 
{v = 0). For improving their effectiveness we 
have used "short memory" principle where 
V = for k < (L/T) or v = k - {L/T) for 
k > {L/T). Besides the "short memory" the 
control quality is influenced by time step T. 

Fractional-order controllers can be realized as soft- 
ware or by passive or active electrical elements [||, |^ . 

6 Example 

We verify the above methods on an example from |^ . 
Assume the system which we can describe by three- 
term [n = 2) differential equation with coefficients 



aa = 0.8, ai = 0.5, aq = l,/?2 2.2, /3i = 0.9, /5o = 0. 
After its approximation with an integer-order sys- 
tem we have 02 ~ 0.7414, ai = 0.2313, oq = 1,/32 = 
2,/3i = l,/3o = 0. The integer-order PD con- 
troller designed with the dominant roots method ap- 
plied to the approximated system has the parameters 
K = 20.5, Td = 2.7343,(5 = 1. By applying the con- 
troller to the original system we do not achieve the 
required quality of the regulation process as with the 
approximated integer-order system. This was con- 
firmed via simulation in the time domain and also 
by checking the stability measure and damping mea- 
sure. 

This proves the inadequacy of approximating the 
fractional-order system with an integer-order sys- 
tem for the purpose of controller design. It is suit- 
able to consider fractional-order systems and also 
controllers, by means of which it is possible to ob- 
tain higher quality regulation and robustness. For 
the fractional-order PD^ controller we then have 
K = 20.5, Td = 5.79,(5 = 0.95 and the required 
quality of regulation is ensured also in the original 
fractional-order system |l^ . 
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Fig. 3: Transient unit - step responses 



On the Fig. 3 is showed the comparison of the unit- 
step responses the integer-order PD controller ap- 
plied on the fractional-order system (thin line) and 
the fractional-order PD^ controller applied on the 
fractional-order system (thick line). 

The dynamical properties of the closed loop with 
fractional-order controlled system and the fractional- 
order controller are better than the dynamical prop- 
erties of the closed loop with the integer-order con- 
troller. The systems with the integer-order controller 
stabilizes slower and has larger surplus oscillations. 
We can see that use of the fractional-order con- 
troller leads to the improvement of the control of the 
fractional-order system. 



7 Conclusion 

The above methods make it possible to design 
fractional-order controllers with given measures of 
stability and damping. The results of previous works 
also show that fractional-order controllers are more 
robust which means they are less sensitive to 
changes of the system parameters and controller pa- 
rameters. This can even lead to qualitatively differ- 
ent dynamical phenomena in control circuits. 
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